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ABSTRACT 

Fingering convection (or thermohaline convection) is a weak yet important 
kind of mixing that occurs in stably-stratified stellar radiation zones in the pres¬ 
ence of an inverse mean-molecular-weight gradient. Brown et ah (2013) recently 
proposed a new model for mixing by fingering convection, which contains no free 
parameter, and was found to fit the results of direct numerical simulations in 
almost all cases. Notably, however, they found that mixing was substantially 
enhanced above their predicted values in the few cases where large-scale grav¬ 
ity waves, followed by thermo-compositional layering, grew spontaneously from 
the fingering convection. This effect is well-known in the oceanographic context, 
and is attributed to the excitation of the so-called “collective instability”. In 
this work, we build on the results of Brown et al. (2013) and of Traxler et ah 
(2011b) to determine the conditions under which the collective instability may 
be expected. We find that it is only relevant in stellar regions which have a rela¬ 
tively large Prandtl number (the ratio of the kinematic viscosity to the thermal 
diffusivity), 0(1CT 3 ) or larger. This implies that the collective instability cannot 
occur in main sequence stars, where the Prandtl number is always much smaller 
than this (except in the outer layers of surface convection zones where fingering 
is irrelevant anyway). It could in principle be excited in regions of high electron 
degeneracy, during He core flash, or in the interiors of white dwarfs. We discuss 
the implications of our hirelings for these objects, both from a theoretical and 
from an observational point of view. 

Subject headings: hydrodynamics - instabilities - stars : interiors - stars : evo¬ 
lution 
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Introduction 


Discovered in the late 1950s in the oceanographic context (Stern 1960), fingering - or 
“thermohaline” - convection was first introduced to the astrophysical community by |Ulrich 


(1972 

) and later revisited by 

Kippenhahn et al. 

(1980) 


The fingering instability occurs 
in stars within radiation zones that have an unstable mean molecular weight gradient (//- 
gradient hereafter). This situation is often found as a result of material accretion onto a star 
by anything from a single or multiple planets ( Vauclair| 2004 Theado & Vauclair 2012), to 
material from a dust-enriched or debris accretion disk (e.g. Deal et al. 2013), or material from 
a more evolved companion (Ulrich||1972 Stancliffe et al.||2007). It also naturally arises in the 


vicinity of the H-burning shell in red giant branch (RGB) stars, as discussed by Charbonnel 


& Zahn (2007) and Denissenkov (2010), and in thin element-rich layers near the surface of 


intermediate-mass stars ( Theado et al. [2009 Zemskova et al.||2014 ). The fingering instability 
initially takes the form of thin tubes, hence the name “finger”, within which the fluid moves 
vertically. The tubes rapidly break down, however, as a result of parasitic shear instabilities 


that develop inbetween them (Radko & Smith 2012 Brown et al. 2013), and the fingering 


instability eventually saturates into a state of homogeneous fingering convection where the 


typical aspect ratio of the eddies is closer to one (Traxler et al. 2011). 


One of the main challenges in the subject of fingering convection is to model the rate at 
which it transports or mixes various quantities of interest, such as heat, angular momentum, 
and chemical species. Until very recently, the only available mixing prescriptions for fingering 


convection were the ones originally proposed by Ulrich (1972) and Kippenhahn et al. (1980). 


Both were essentially based on dimensional arguments and ad-hoc modeling efforts, that 
could not, for lack of experimental data, be directly verified. Furthermore, the two estimates 
for the turbulent diffusion coefficient differed by several orders of magnitude. 

Thanks to advances in supercomputing, however, much progress has recently been made 
to measure the rate of mixing in fingering convection in numerical experiments, and to de¬ 


velop theoretical models for the collected data. Traxler et al. (2011) first proposed a sim¬ 


ple empirical formula for the turbulent diffusion coefficient for both heat and composition, 
which was fitted to a substantial number of three-dimensional numerical experiments run at 
various input parameters (notably, different values of viscosity, thermal diffusivity, compo¬ 
sitional diffusivity and background stratifications of temperature and composition). Brown 


et al. (2013) later added further simulations and showed that the model of Traxler et al. 


(2011), while perfectly adequate for more strongly stratified systems (ie. systems with a 


weaker inverse ^-gradient, or systems with stronger thermal stratification), underestimates 
the mixing efficiency for systems that are very close to being unstable to standard over¬ 
turning convection. They proposed an alternative model, based this time on first principles, 
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which fits the numerical data remarkably well for most available simulations (in the limit 
where the compositional diffusivity is smaller than the kinematic viscosity, which is always 
the case in stellar conditions), and is simple enough to be used in real time in any stellar 
evolution code. 

Despite this significant progress, much remains to be done to fully characterize the 
properties of fingering convection in astrophysics. In the related oceanographic context, it is 
well-known that small-scale homogeneous fingering is itself unstable to secondary large-scale 
instabilities of various forms including the 7 -instability, which leads to the formation ther¬ 
mohaline staircases (Radko|[2003), the collective instability, which excites large-scale gravity 


waves (Stern 1969 Stern et ah 2001), and the intrusive instability, which also gives rise to 


layering in the presence of horizontal thermo-compositional gradients (Holyer 1983 Walsh 


& Ruddick 1995). In all three cases, the presence of larger-scale dynamics can significantly 


enhance the rate of turbulent mixing above that of homogeneous fingering convection (see 
Brown et al. 2013, for instance), so it is crucial to understand why and when these secondary 
instabilities arise. 


As shown by Traxler et al. ( 2011 b), all three large-scale modes of instability can be 
studied within the same unifying framework, which consists in deriving so-called “mean- 
held” equations for the large-scale fields while treating the small-scale fingering convection 
only through its contribution to the turbulent transport of heat, composition and momentum. 
However, their work was specifically applied to oceanographic conditions; whether similar 
secondary instabilities also occur in stellar conditions remains the subject of undergoing 
research. It has recently been shown that the 7 -instability cannot operate in fingering systems 


in stellar interiors ( 

Traxler et al. 2011 

Brown et al. 2013 

), but that intrusive instabilities are 

ubiquitous in the presence of horizontal gradients ( 

Medrano et al. 

2014 

). In this paper, we 


focus on the third form of large-scale instability, namely the collective instability. 


The collective instability was first studied by Stern (1969) in the oceanographic context 


(see also Stern et al. 2001). As mentioned above, it is a mechanism by which small-scale 


fingering convection naturally excites large-scale internal gravity waves. By “large” we imply 
waves whose typical wavelength is at least one order of magnitude greater than that of the 
basic fingering instability. The physical mechanism behind the collective instability is easy 
to understand, given some basic knowledge about the development of the standard double- 


diffusive instabilities, which we now briefly summarize (see the review by Garaud (2013) for 
more detail). First, recall that double-diffusive instabilities in general develop in fluids whose 
overall density is stably stratified, but depends on two different quantities that ( 1 ) diffuse at 
different rates and (2) have opposing contributions to the density stratification. The fingering 
instability arises whenever the faster-diffusing component (typically, temperature) is stably 
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stratified while the slower-diffusing (typically, composition) component is unstably stratified. 
The oscillatory instability (which is related to semiconvection), by contrast, develops when 
the faster-diffusing component is unstably stratified, while the slower-diffusing component is 


stably stratified (Kato 1966 Baines & Gill 1969). This occurs in thermally-unstable regions 


that are stabilized by a strong ^-gradient. The oscillatory instability takes the form of 
overstable gravity waves, i.e. oscillatory motions whose amplitude increases exponentially 
with time. 


As clarified by Radko (2013), the collective instability can simply be viewed as a large- 


scale manifestation of the oscillatory instability, in a system which (on the small scales) 
is unstable to fingering convection and where the microscopic thermal and compositional 
diffusivities are therefore augmented by the induced turbulent mixing. Indeed, the key is to 
note that fingering convection transports composition much more rapidly than it transports 
heat, so the turbulent diffusivity of the stably stratified field (composition) is much larger 
than the turbulent diffusivity of the unstably-stratified field (temperature). This situation, 
as discussed above, gives rise to an oscillatory instability, whose growth rate depends on 
the fingering fluxes. Furthermore, since the turbulent diffusivities are larger than their 
microscopic counterparts, the typical lengthscale of the collective modes is also significantly 
larger than the basic fingering scale. This has advantages and disadvantages for modeling 
the instability. On the one hand, the separation of scale between the collective instability 
and the fingering instability allows for the use of mean field theory, as shown by |Traxler et al. 


(2011b), a particularly convenient tool in this problem. On the other hand, it also means 


that the collective instability can only be studied numerically in very large computational 


domains, which has only recently become feasible. Stellmach et al. (2011) were the first to 


demonstrate the natural emergence of the collective instability in numerical simulations of 
fingering convection at parameter values similar to those expected in the ocean. They found 


that the theory developed by Traxler et al. (2011b) adequately predicts the observed growth 


rate of the collective modes in that case. 


Denissenkov & Merryfield (2011) were the first to study the collective instability in 
the astrophysical context, with particular application to RGB stars. Using the formalism 


of Traxler et al. (2011b) together with estimates of the fingering fluxes from 2D numerical 


simulations, they found that collective modes are not excited. On the other hand, Brown 


et al. (2013) saw the emergence of oscillatory features in their 3D direct numerical simulations 


of fingering convection. This apparent discrepancy can be explained by the difference in 


the parameters used: the results of Denissenkov & Merryfield (2011) were obtained for 
realistic stellar values of the Prandtl number and diffusivity ratio (Pr = u /kt ~ 10 -6 and 
r = k^kt r^j 10 ' respectively), while the simulations of Brown et al. (2013) were run, 
for numerical feasibility, at Pr ~ r ~ 0.01 instead. The difference between the findings 
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of Denissenkov & Merryfield (2011) and Brown et al. (2013) raises an interesting question, 
however: where in parameter space does the collective instability lie, and are there any stars 
for which it may be relevant? 


The possibility of driving large-scale gravity waves in fingering convection is indeed 
quite exciting for two reasons. First, and as discussed above, these secondary instabilities 
are usually accompanied by enhanced transport, and could therefore help solve some of 
the outstanding problems in stellar evolution associated with “missing mixing”. Second, 
the gravity waves themselves may be directly detectable via asteroseismology, should they 
have sufficiently large spatial structures and amplitudes. Our goal, therefore, is twofold: 
to perform a comprehensive study of the collective instability and determine under which 
conditions (i.e. for which values of the microscopic diffusivities, background thermal and 
compositional stratification) the latter is expected, and then to search across the ffR diagram 
for stellar objects where such conditions may be realized. 


In what follows, we study the collective instability in the more general astrophysical 
context (i.e. without limiting our study to parameters relevant for RGB stars) in Section [2j 
We also improve on the study by Denissenkov & Merryfield (2011) by using the more recent 
fingering ffux laws of Brown et al. (2013), and by comparing our results with large-scale 3D 
numerical simulations of the collective instability (see Section [3j. As we shall demonstrate 
the latter is only present for values of the Prandtl number and of the diffusivity ratio down to 
about 10 -3 , in regions of parameter space that are very close to being overturning-unstable. 
This finding is therefore consistent both with the results of Denissenkov & Merryfield (2011) 
on the absence of the collective instability in RGB stars and on its presence in the numerical 
simulations of Brown et al. (2013). We then look in Section [4] across the HR diagram and 
attempt to identify stars which have “large” Prandtl number and diffusivity ratio. Our 
findings are discussed in Section |5j 
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2. The collective instability 


To model the collective instability, we begin with the standard (non-dimensional) equa¬ 
tions of fluid dynamics applied to the problem of fingering convection (Radko 2003): 


1 f d\i 

idw +u ' Vu 


= -Vp+ (T - p)e z + V 2 u, 
Vu = 0, 


dT 


V + U-VT + W = V 2 T, 
at 

dg „ 1 ,, 

vrr + u • V/i + —w = rV 2 fi, 
at Ro 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


where u = (u, v, w) is the velocity field, and p, T and /j are the pressure, temperature and 
mean molecular weight perturbations away from the chosen background. Assuming that 
the vertical size of the domain under consideration is small compared with the respective 
local scale heights, the background temperature and mean molecular weight profiles are 
approximated by linear functions of the vertical coordinate z, with constant (dimensional) 
gradients T 0z and p 0 z respectively. To arrive at the non-dimensional equations 0-0 , we 
have used the following system of units. The unit length is 


[l]=d 


f __KtO__\ 1/4 
\ga(T 0z -T$)J 



( 5 ) 


where Kt is the thermal diffusivity, v is the viscosity, g is gravity, a is the coefficient of 
thermal expansion, is the local adiabatic temperature gradient, and Nt is the local 
buoyancy frequency based on the thermal stratification only. All of these quantities are also 
assumed to be constant. The unit time is [t] = cP/kt, and the unit velocity is [u] = Kr/dt. 
The unit temperature is [ T] = d(T 0z — Tfif) and finally, the unit mean molecular weight 
is [/i] = ad(T 0z — where f3 is the coefficient of “compositional contraction”. The 

coefficients a and /3 are thermodynamic derivatives of the equation of state, and given by 




( 6 ) 


where po is the mean density of the region considered. The three standard non-dimensional 
parameters of fingering convection thus emerge: the Prandtl number Pr = the diffu¬ 
sivity ratio r = and the density ratio Ro = a(T 0z — T^)/j3po z - A system is linearly 

unstable to fingering convection when Rq is in the range [l,r _1 ] (Stern 1960 Baines & Gill 


1969), so we restrict our following analysis to that range only. 
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We now filter (jl])-(|4]) to extract the dynamics of scales that are significantly larger than 
the basic finger scale (which is of the order of 10d, typically). This leads to 

-Wp+(T-Jl)e z + V 2 u, (7) 

= 0 , ( 8 ) 

= V 2 T, (9) 

= t\7% ( 10 ) 

where the overbar denotes a local spatial average over the small fingering scales, that we 
assume commutes with both spatial and temporal derivatives. Defining u' = u — u and 
similarly for T' and //, averages of the nonlinear terms in the original equations give rise to 
the Reynolds stress tensor R = u'u', the temperature flux F T = n'T' and the compositional 
flux Fj, = u'/i'. This system of equations is not closed unless these three quadratic correlation 
terms are known functions of the large-scale variables u, T and ~jJ and of the model parameters 
Pr, r and Rq. 


1 f du 

Pr \dt + n ' Vn + V ' R 

Vu 

df _ - __ 

——h u • VT + V • F t + w 
at 

djl _ __ 1 _ 

_ +u . V ^ + V.F„+- ro 


We now follow the procedure proposed by Traxler et al. (2011b) to model the collective 
instability. The steps are outlined here for completeness; the reader is referred to the original 
paper for more pedagogical detail if needed. Traxler et al. (2011b) proposed the following 
“closure” model: (a) the Reynolds stress is negligible compared with other terms in the 
momentum equation, (b) the turbulent fingering fluxes are principally vertical, so Ft = FVe 2 
and similarly for F /t = F^e z . Condition (a) has been verified a posteriori via numerical 
simulations. Condition (b) on the other hand is only marginally justified as numerical 
simulations show that the horizontal components of the fluxes can be up to a third of the 
vertical component. Nevertheless, we shall make that assumption for simplicity. We now 
define the thermal Nusselt number Nu^ and the turbulent flux ratio 7 as 


F7 1 — (1 — Nuj 1 ) 



( 11 ) 


In the final assumption of the closure model, (c), we assume that Nut and 7 are functions 
of T and ~jJ only through the local density ratio 


R 


Rq 


i + S 
1 + -R0S' 


( 12 ) 


Nu t and 7 of course also depend on Pr and r, so (c) implies Nu^ = Nu t(R] Pp t) and 
7 = 7 ( 7 ?; Pr, r). Whether assumption (c) is valid or not will be verified a posteriori, if the 
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collective instability model correctly predicts the growth rate of large-scale gravity waves 
observed in simulations. We do know, however, that (c) is valid in the absence of large-scale 
temperature and compositional perturbations, since in that case the fluxes can only depend 
on the input non-dimensional parameters of the system, i?o, Pr and r. 


Substituting GD and ( Jl2| ) into (|7|)-(fl0]), we immediately see that there is a trivial 
solution to the resulting equations: one where all the large-scale fields are zero, namely 
u = T = ~jl = p = 0 , where the local density ratio is equal to the background density ratio, 
R = Ro, and where Ft and F tJ are constant and equal to their corresponding values in the 
absence of large-scale perturbations, namely 


Ft = 1 — Nut(-Ro; Pr, r) = 1 — Nuo , 
_ F t _ F t 

“ 7(-Ro;Pr,T) - 70 ’ 


(13) 


which defines Nuo and 70 . This solution is merely the state of homogeneous fingering con¬ 
vection one would expect in the absence of large-scale perturbations. 


We now consider that T, jj and u are large-scale but small-amplitude perturbations 
on this homogeneously turbulent state, and linearize the system ([T])-( 10 ) around the latter. 
Linearizing (12) around Rq, the local density ratio becomes 


R 


Rq 



-R 


<¥\ 

°dz ) 


(14) 


and we define the coefficients A\ = i?o (^7 l /dR)R=R 0 and A 2 = Ro(dNi\T/dR) r=r 0 such 
that 

T\R-,Pi,t) c To" + (R- fio) - TcT 1 +-4. (f- «o|) . 

Nu t(R] Pr, r) ~ Nu 0 + {R- R 0 ) — Nu o + A 2 (jp ~ R o^j ■ ( 15 ) 

Assuming that all quantities can be expressed as linear combinations of normal modes of the 
kind q = qexp(ilx + imy + ikz + At), we then obtain a cubic equation for the growth rate 
A of the large-scale modes, i.e. 


A 3 -|- g^A^ T cl\ A -|- oq — 0 , 


(16) 
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where 


a 2 = |k| 2 (l + Pr + r) + k 
ai = jk| 4 (rPr + r + Pr) + fc 2 |k |J 


(1 — A!i?o)(NUo — 1) + A 2 ( 1- 

7o 


(17a) 


(r + Pr)(A 2 + Nu 0 - 1) - A 2 (l + Pr)^ 

7o 


-Aii? 0 (l + Pr)(Nu 0 - 1) 


- fc 4 Aii?o(Nuo - l ) 2 + Pr 


2 L 2 f l 

2 , - ( l - 

kl 2 V R 


(17b) 


ao — 


|k| 6 rPr + fc 2 |k| 4 Pr 
L 2 


(r - Aii? 0 )(Nu 0 - 1) + A 2 ( r - — 

7o 


o 

— /c 4 |k| 2 Pri?o^4i(Nu 0 — 1)" 


+Pr- 


kP 


T- — ) + k 2 A\(l — i? 0 )(Nu 0 — 1) 


-k 2 [A 2 (1 - Ro) + Nuo - 1] ( -- 

R) 7o 


(17c) 


where |k | 2 = l 2 + m 2 + k 2 . This is the same cubic as in Traxler et ah (2011b), assuming that 


there are no background horizontal gradients of temperature and composition. Note that 
in this limit, there is a perfect symmetry between the x and y directions, and the system 
dynamics only know about the total horizontal wavenumber L = y/l 2 + m 2 . Also note that 


the cubic (16) has 3 solutions A for each input k. In what follows, when referring to “the 


growth rate of the collective mode with wave vector k” we always implicitly take the root 
which has the largest real part. 

The only remaining problem is to determine the constants Nuo, 7 o, A\ and A 2 , which are 
needed to compute do, «i and a 2 . They are determined from the properties of homogeneous 
fingering convection, and are functions only of the density ratio Rq and of Pr and r. While 


Traxler et al. (2011b) used values of these coefficients appropriate for fingering convection in 
an oceanographic context (Pr ~ 7 and r ~ 0.01), in stellar interiors Pr and r are typically 
much smaller than one, so their values cannot be used here. However, Nuq and 70 can be 


directly estimated using the model proposed by Brown et al. (2013) for mixing by fingering 
convection: 


Nu t (R', Pr, t) = 1 + 49 


X 2 

Agm 


7 (i?; Pr, r) = R 


^fgm^fgm “I" 4gm) 
' fgm "b T ^fgm 


Afgm + l{ gm 


(18) 


for R = R 0 , where Af gm and lf gm are functions of R, Pr and r and are the growth rates and 
horizontal wavenumbers respectively of the fastest growing modes of the basic fingering insta¬ 
bility at these input parameters (recalling that, for the basic instability, the fastest-growing 
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modes have a zero vertical wavenumber). They can be found by solving simultaneously a 
cubic and a quadratic equation: 

Af gin + ^ 2 ^fg m + b\ Afg m + &o = 0, where (19) 

i>2 = ^f gm (l + Pr + r), 

b i = 4m( rPr + Pr + t) + Pr ^1 - -0 , 
b 0 = /| m rPr + /| m Pr (r - , 


c 2 Af gm + CiA fgm + c 0 = 0, where 
c 2 = 1 + Pr + r, 
ci = 2/| m (rPr + r + Pr), 

c 0 = 3Zf gm rPr + Pr i T ~ ^ 


( 20 ) 


Note that (19) is simply obtained from (16) with q 0 = Nuq = A 1 = A 2 = 0 (which recovers 


the dispersion relation for the basic fingering instability), R 0 = R and k = 0 , while ( 20 ) 


is then derived from (19) by maximizing the growth rate A over all possible values of the 


horizontal wavenumber. The coefficients Ai and A 2 are then obtained by differentiating the 
functions Nu(r(-R;Pr, r) and 7 _ 1 (i?;Pr,r) numerically, as 


A\ — R( 


A 2 — Rt 


7 1 (Rp + e; Pr, r) — 7 1 (R 0 -,Pt,t) 
e 

Nu T (i ? 0 + e; Pr, t) - Nu t(Rq) Pr, r) 


( 21 ) 


In practice, we take e = 10 3 Rq in what follows. 

Figure [I] shows the real part of the growth rate Re (A) of each mode with spatial structure 
given by a wave vector of horizontal component L and vertical component k, for various values 
of the Prandtl number (assuming r = Pr), and for various values of Rq- Modes with negative 
Re( A) are shown in white. Note that for ease of comparison of the various diagrams, we use 


the reduced density ratio introduced by Traxler et al. (2011), namely 


r = 


Ro — 1 


- 1 - 1 ’ 


r 


( 22 ) 


which is 0 when R 0 = 1, i.e. when the system is marginally unstable to standard overturning 
convection, and is 1 when R 0 = t^ 1 , ie. when the system is marginally stable. Hence cases 
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with r = 10 4 are close to being overturning-unstable, while cases with r = 0.1 are only 
weakly unstable to fingering. 


We immediately note that the growth rate is independent of the sign of k, which is 
expected since the cubic governing A only depends on k 2 . The same is true for L , so the 
figure only shows results for positive values of the latter. As found by Traxler et al. (2011b), 
the instability diagram thus drawn takes the shape of a flower. The bulb of the flower 
centered around L ~ 1 and k — 0 corresponds to the basic fingering modes, while the leaves, 
when they exist, are the large-scale gravity-wave modes excited by the collective instability. 
By contrast with the fingering modes, the collective modes have a complex growth rate, 
whose imaginary part is the oscillation frequency of the wave. They typically have low L 
and low (but non-zero) k wave numbers. Figure [I] suggests that they only exist for the 
smaller values of Rq and the larger values of the Prandtl number (and r). 

To see this more clearly, we show in Figure [2] the growth rate (real and imaginary parts) 
and wavevector of the fastest-growing collective mode as a function of Pr, r and of the 
reduced background density ratio r. These are found by maximizing the real part of A over 
all modes of wavenumber L and k with non-zero imaginary part. We see that collective 
modes exist for Pr = r = 0.1, 0.01, and 0.001, but not for Pr, r < 10~ 4 . We also see that for 
Pr, r > 10 -3 , these modes only exist for very low values of r (equivalently, for density ratios 
very close to one), in other words, for systems whose background stratification is already 
fairly close to being Ledoux-unstable. Note that this result had already been obtained by 


Brown et al. (2013) using the Stern number as a approximate diagnostic of the presence of 


the collective instability. 


The reason why the collective instability is limited to the higher values of Pr and r, and 
the lower density ratios, is quite simple: this also corresponds to the region of parameter 
space where the magnitudes of the fingering fluxes (of temperature and composition) are 
large enough to have a significant effect on the overall heat and compositional transport 
when compared with pure diffusion. Since the physical mechanism behind the collective 
instability discussed in Section 1 crucially relies on the turbulent fluxes, it is not surprising 
to find that the instability disappears when these fluxes are too small. 


Wherever collective modes exist, we find that the real part of A for the fastest-growing 
ones is of the order of Pr, which in dimensional terms corresponds to Pr k^/c/ 2 = Pr 
For Pr ~ 10~ 2 , this is of the order of 0.1A^. The imaginary part of A is of the same order, 
so we expect each of these modes to oscillate with a frequency commensurate with their 
exponential growth rate. Finally, both horizontal and vertical wavenumbers are of the order 
of a few times 1CT 2 , so we expect the fastest-growing collective modes wavelengths to be of 
the order of a few hundred times d or in other words, a few tens of times the width of the 
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fastest-growing fingering mode. What this corresponds to in real dimensional terms depends 
on the actual dimensional values of the thermal and compositional diffusivities and of the 
local buoyancy frequency, and will therefore vary from star to star (see Section [5] for an 
example). 


3. Comparison with numerical simulations 


In order to test the validity of the collective instability theory described above, we now 
turn to Direct Numerical Simulations (DNS) of fingering convection, and use a computational 
domain that is sufficiently large to contain the large-scale gravity waves we expect will 


develop. We use the same code as Traxler et al. (2011b), which solves the set of equations 


(JTJ-Q in a triply periodic domain. We present here the results of a simulation at Pr = 0.1 
and t = 0.03, which are the smallest values of the Prandtl number and of the diffusivity 
ratio we can realistically achieve while at the same time using a computational domain large 
enough to contain at least a few wavelengths of the fastest-growing collective mode. We use 
a density ratio of Rq = 1.1, for which the growth of that collective mode is relatively fast. 
This is the same set of parameters as one of the simulations of Brown et al. (2013) for which 
wave-like behavior (and eventual layer formation) was indeed observed. 

At these parameters, the non-dimensional growth rate and horizontal wavenumber Af gm 
and /f gm of the fastest-growing fingering mode are: 


A fgm ~ 0.2 and / fgm ~ 0.75, 


(23) 


while the growth rate, horizontal wavenumber and vertical wavenumber of the fastest-growing 
collective mode are 


A coll ~ 0.024 ± O.IOz , Leon ~ 0.060 and k con ~ 0.048 . 


(24) 


As a result, we need a domain whose vertical extent is at least 27r/0.048 ~ 130d, with a 
similar horizontal extent, to contain the fastest-growing mode of the collective instability. 


For this reason, we cannot use the aforementioned simulation of Brown et al. (2013), which 


was performed in a domain of size (100d) 3 . In what follows, we use a domain of size ( L x = 
400 d, L y = 20 d,L z = 400d), which is very thin in the y-direction to save on computational 
time. It is nevertheless thick enough to contain at least two wavelengths of the fastest- 
growing fingering mode (which is about 8d at these parameters). This was found to be 
an excellent trade-off that captures the three-dimensional dynamics of the basic fingering 
instability adequately while being quasi-two-dimensional for the collective instability. The 
effective resolution in this run is 384 meshpoints per lOOd. 
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The simulation is initialized with small random fluctuations in the temperature field, 
which rapidly grow according to the dynamics of the fingering instability. This is illustrated 
in Figure [3j which shows the total kinetic energy in the domain as a function of time. As 
described by Brown et al. (2013), this initial linear instability then saturates as a result 
of the shear that develops between the fingers. This saturation can be seen in the first 
plateau observed shortly after t = 90. Figure [4|a shows a snapshot of the compositional 
perturbations at t — 100, revealing a fairly homogeneous fingering field. Shortly thereafter, 
various large-scale gravity waves excited by the collective instability begin to grow. Figure 
i> shows the same simulation as in Figure [4^,, but at a later time when collective modes 
are clearly visible. Figure [3] shows the kinetic energy in some of the largest-scale collective 
modes, and compares the observed growth rates with the corresponding ones predicted by 
the theory outlined in Section [2} We find that our predictions are typically quite good, which 
confirms the assumptions made in Section [2] 


The collective instability appears to saturate when the kinetic energy in individual 
modes reaches a small fraction (a few percent) of the total kinetic energy. However, given 
that each mode is part of a mode family^ whose members all grow at the same rate and 
saturate at the same amplitude, saturation occurs in fact when there is roughly equipartition 
between the energy in all the fingering modes and the energy in all the collective modes. 


As seen in Figure [5j the growth of the waves causes a gradual increase in the global 
mixing rate. This can be quantified through the compositional Nusselt number, which is 
the ratio of the effective turbulent diffusivity D^ to the microscopic diffusivity of the mean 
molecular weight: 


Nu m = ^ 

K/j, 


(25) 


Interestingly, we see that Nu„ continues to increase even when the wave amplitudes no 


longer do. This signals the onset of layered convection and confirms the results of Brown 


et al. (2013), who also found that their simulation transitioned into layers for the same 
set of parameters. Layered convection is characterized by a horizontally-averaged density 
profile that takes the form of a staircase, with region of near-uniform density (the convective 
layers) separated by strongly stable interfaces. We see in Figure [6] that two layers form, 
with a spacing similar to the wavelength of the fastest-growing collective mode. This is 
consistent with the results of Traxler et al. (2011b) in the oceanographic context, who also 
found that emerging thermo-compositional layers have the same wavelength as that of the 
fastest-growing collective modes in their simulations. By contrast with the high-Prandtl 
nature of oceanic fingering, the layers are not very robust, and are constantly pierced by 


1 a mode family being all modes of the kind ( ±k x , 0, ±fc 2 ) 
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strong updrafts and downdrafts that ultimately control the chemical transport rate in this 
system. By the end of the simulation (around t = 1000), we begin to see hints that the 
two layers may be merging as a result of these strong mixing events. In order to study 
the overall dynamics and general transport properties of layered convection in this regime, 
however, numerical simulations in domains of more even aspect ratio are needed. Indeed, 
at this point the quasi-two-dimensional nature of the simulation presented here becomes a 
hindrance as the convective eddies are severely restricted by the aspect ratio of the domain. 
Given the resolution required, this run would effectively have to be a (1500) 3 run, which is 
beyond the scope of this paper. 


It is worth noting that the mechanism by which thermo-compositional layers form in 
our simulation remains to be determined. Indeed, the standard layering transition caused 


by the 7 -instability cannot occur at this parameter regime, as demonstrated by Traxler 


et al. (2011). An alternative possibility is that these layers are forming as the result of the 


nonlinear development of the collective modes, as originally discussed by Stern et al. (2001) 
in the oceanographic context. This scenario is much more difficult to confirm, since it can 
only be studied using DNS. Hence, whether thermo-compositional layering is an inevitable 
consequence of the collective instability, or whether it only occurs in a subset of parameter 
space remains to be determined. 


4. Application throughout the HR diagram 

We now look throughout the HR diagram to see whether any stars harbor regions of 
relatively high Prandtl number and diffusivity ratio. We focus on cases that are not fully 
convective, thereby ignoring very low-mass stars and stars on the Hayashi track. We split 
our study between stars on the main sequence, and post main sequence objects. 


4.1. Main sequence stars 


The presence of inverse yU-gradients is fairly rare in main sequence stars. This is because 
nuclear reactions, which tend to increase the mean molecular weight, are more efficient at 
higher temperatures, hence closer to the core. The gravitational settling of He with respect 
to H during the main sequence also usually contributes to the formation of a strong stable 
/i-gradient. There are, however, three notable exceptions to this general statement. The first 
is the possibility of 3 He burning, which is the only common reaction that decreases the mean 
molecular weight of the material. As first discussed by Ulrich (1971) and Ulrich (|I972), this 
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effect can drive fingering convection, and can be important in RGB stars for instance (see 
section 4.2). The second exception is through the formation of thin element-accumulation 
layers near the stellar photosphere through the combination of radiative levitation and grav- 


itational settling ( 

Michaud 

1970 

Charpinet et al. 

1997; 

Richard et al. 

2001 

). These layers 

can also be unstable to fingering convection, as first discussed by 

Theado et al. 

(200fF 

) and 


more 


Finally, fingering convection can also be triggered by the accretion of higher-metallicity 


material on the surface, either via planetary infall ( 

Vauclair 

2004; 

Garaud 

2011 

) or from a 

more evolved companion (jUlrich 

1972 

Stancliffe et al. 2007 

). For stars roughly below 1.5M 0 , 


the accreted material is first rapidly mixed within their outer convection zone (which may 
deepen slightly as a result of the added metallicity). This creates an inverse //-gradient at 
the bottom of the convection zone, which can be unstable to fingering convection under 
the right circumstances. For stars above 1.5M 0 , which normally have a radiative envelope, 
the accreted material may at first trigger the formation of a shallow outer convective layer. 
The latter rapidly mixes the added material, then disappears leaving behind a strong /i- 
gradient, which can drive a fingering instability. The possibility of subsequently triggering 
the collective instability in the process has however never been investigated. 


We now look at stars in the mass range O.7M 0 -3OM 0 , and determine their Prandtl 
number and diffusivity ratio. Since we are only interested in an order of magnitude estimate, 
we approximate the thermal diffusivity by its radiative contribution: 


HvT ^rad 


16crT 3 
3 np 2 c p ’ 


(26) 


where a is the Stefan-Boltzmann constant, T is the temperature, k is the opacity, p is the 
density, and c p is the specific heat at constant pressure. We estimate the viscosity from the 
sum of its radiative and collisional contributions: 


16crT 4 0A06m]i 2 (k B T) 5 / 2 

" _ "“ d + ~ + eV(hiA„ H )p 


(27) 


where c is the speed of light, k B is the Boltzmann constant, e is the electron charge, m B 
is the mass of the hydrogen atom, and In Ay is the Coulomb logarithm for the collisions 
between element i and element j: 


In Ay 


-19.26 


- In p + - In T — - In 
2 2 2 


X+ 3 


In (ZiZj) , 


(28) 


where Zi is the charge of element i, and Zj is the charge of element j. The expression for 
zz co n is from Spitzer (1965). The function (7(111 A) is a correction to the Coulomb logarithm 
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suggested by Michaud & Proffitt (1993), and is given by 


<7(111 A) 


1 

L2 


In (exp(1.2 In A) + 1) . 


(29) 


For simplicity, when estimating the viscosity we use a fully ionized pure hydrogen gas, in 
which case Z, = Zj = 1 and X = 1. The correction to the Coulomb logarithm for a H-He 
mixture is of the order of a few percent, and is neglected here. These estimates for kt and 
v are valid as long as the stellar material is non-degenerate. 


Finally, we use the diffusion coefficient given by Michaud & Proffitt (1993) 


15 /2 m H ( k B Tf / 2 3 3 + X 

^ ~ 16pIn A HHe V "ffir e 5 (1 + X)(3 + 5X)(0.7 + 0.3X) ’ 


(30) 


for the diffusion of He in a H gas (or vice versa). This coefficient overestimates the one for 
the diffusion of heavier elements in a H-He mixture^ but the difference does not affect our 
estimate of where the collective instability may occur. 


The Prandtl number Pr = v/K t-, and the diffusivity ratio r = Kt thus calculated are 
shown in Figure [7] as a function of radius for stars of different masses. The stellar models 
used to calculate these ratios are obtained using the MESA stellar evolution cod^] (Paxton 


et al. 2011, 2013) to evolve each star from zero-age main sequence until their central mass 


fraction of hydrogen drops below 0.4. In each case, the star has initial solar composition, 
and the only mixing considered was due to standard convection. Note that while we have 
picked a particular stellar evolutionary stage at which to present this data, it is a good 
representation of both Pr and r for these stars across the entire main sequence. We see 
that the diffusivity ratio is always smaller than the Prandtl number, by a factor of a few in 
the lower-mass stars (where v is dominated by collisions), and by a much larger factor in 
high-mass stars (where it is dominated by radiation). Generally-speaking, we also see that 
the only regions which contain material with a large Prandtl number are near the surface 
layers of stars of mass lower than 1.3M 0 , which corresponds precisely to the outer parts 
of their convective regions which cannot support fingering. We therefore conclude that the 


2 A very rough order-of-magnitude estimate of the diffusion coefficient of a trace species in a H-He mixture 
can be made by dividing k^ given in equation (30) by Z 2 , where Z is the atomic charge of the species 
considered (Michaud & Proffitt 1993 ). Using equation ( 30 1 instead therefore overestimates the true diffusivity 
by a factor of a few to a few hundreds depending on Z. As we show later, k^ as calculated using ( 30 ) is 
already too small to give rise to the collective instability. Accounting for the Z correction can only reinforce 
our conclusion rather than change it. 
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collective instability is unlikely to play any role in main sequence stars of any mass or age, 
and that the model proposed by Brown et al. (2013) remains a good estimate for the rate of 
mixing by fingering convection wherever it occurs. 


4.2. RGB stars 


Fingering convection is more common after the main sequence turn-off, during off-center 
(shell-) burning phases. This is the case for instance in RGB stars, where 3 He-burning on 
the outer edge of the H-burning shell creates a small, albeit significant, inverse /i-gradient 


(Eggleton et al. 2006). The effect of fingering convection in that region has been explored 


in depth by Charbonnel & Zahn (2007) and Denissenkov (2010) for instance. Later on, 


fingering convection is also expected during the core He flash, as first discussed by [Ulr ich 
(1972). There, the off-center burning of He generates higher-// material in a thin shell 


surrounding the core, and the lower part of this shell can become fingering-unstable. 

The core of RGB stars is mostly composed of pure fully ionized He, where the electrons 
are degenerate (or at least partially degenerate) while the nuclei remain non-degenerate. In 
this case, the thermal conductivity and shear viscosity are dominated by electron conduction, 
while the compositional diffusivity remains dominated by collisional processes between the 
nuclei themselves. We use the formulae for the shear viscosity r] and thermal conductivity 


kx of strongly degenerate electrons given by Hubbard (1966), namely 


V = 


8 H 5 


5/3 


135 TT 2 mle 4 Z 2 (^3 Amn ) ~ F ’ 

( 27 t h) 3 kg 


kx — 


16 mie^Am 


-Gx(kf)pT . 


(31) 


H 


where h is the Planck constant, A is the average atomic number of the nuclei, m e is the elec¬ 
tron mass, and Kf is the dimensionless Fermi wavenumber kf = (97TZ/4) 1 / 3 . The functions 
Hr and Gr are 


Gr(x) — 


and Hr(x) = 


G r (x) 


where, finally, 


ln(l + f) v ~' 2-2G r (*) + g’ 

Z 2 e 2 / Anrn \ 1/3 


r = 


/ 47T/} 


(32) 


(33) 


ksT \3ArriH, 

These expressions are only valid in a strongly degenerate limit, in which the degeneracy 
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parameter C 2 is small, with 

The thermal diffusivity and kinematic viscosity due to electron conduction in this highly 
degenerate limit are then given by v e = p/p and K e = kx/pc pi as usual. For the purpose of 
the following calculation we use A = 4 and Z — 2 since the core is mostly composed of He. 

For simplicity, because we are again only interested in an order-of-magnitude estimate 
of the Prandtl number and of the diffusivity ratio, we compute the total thermal diffusivity 
and kinematic viscosity by adding the degenerate electron contribution to the usual non¬ 
degenerate terms, as 


K-T — ^rad + j 
v = u rad + Z/coii + V e . 


(35) 


The non-degenerate nuclei contributions typically dominate in the star’s envelope, while the 
electron conduction contribution dominates in the core. Simply adding the two contributions 
is adequate either in the highly-degenerate or in the non-degenerate limits, the relative error 
committed being largest (and of order unity) in the intermediate region of weak degeneracy 
where C 2 ~ 1. Finally, since the nuclei are not degenerate, we use the same formula for k p 
as the one given in equation (30), with liiAnHe replaced by lnAnec- 


The profiles of the Prandtl number Pr and diffusivity ratio r in a 1 M 0 stellar model 
are shown in Figures [8] and [9] at several points on the RGB. Figure [8] shows a profile taken 
just before the luminosity bump, which corresponds to the onset of fingering convection in 
low-mass stars (Charbonnel & Zahn 2007). Figure [9] shows a number of profiles taken during 
the first core He flash at the tip of the RGB. These models were obtained using the MESA 
stellar evolution cod^jto evolve a 1M 0 nonrotating star with metallicity Z = 0.02 starting 
from the pre-main sequence, including atomic diffusion and gravitational settling. 


In both Figures, we see that the Prandtl number and the diffusivity ratio are much 
larger in the core where microscopic transport is dominated by electron conduction. By 
contrast they are quite small in the non-degenerate envelopes, as it is the case for main 
sequence stars. This suggests that the only potential place where the collective instability 
could be at play is in regions of strong electron degeneracy. Having a region of relatively high 
Prandtl number and diffusivity ratio is not enough, however, to guarantee the existence of 
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the collective instability. This region must also be the subject of active fingering convection, 
and the density ratio must be sufficiently low as discussed in Section [2} 


For stars ascending the RGB, the inverse //-gradient begins on the outskirt of the H- 
burning shell, distinctively outside the core. For the 1M 0 star shown in Figure [8j this 
corresponds to the region between m/M 0 = 0.24 and the bottom of the convection zone. We 
see that the Prandtl number is lower than 10~ 4 . Furthermore, the edge of the core itself is the 
seat of a very strong stabilizing /i- gradient, preventing any of the turbulence induced by the 
fingering instability from penetrating into the degenerate region. We therefore confirm the 
results of Denissenkov & Merryfield (2011), namely that the collective instability is unlikely 
to play a role in these stars. 


The case of stars undergoing core He flash is a priori more promising: the initial flash 
takes place well-within the region of strong electron degeneracy where the Prandtl number is 
roughly 10 -3 , a value for which the collective instability is indeed active. The off-center He 
flash convection zone converts a few percent of its 4 He by mass into 12 C via the triple-alpha 
process, leaving an inverse //-gradient at its bottom boundary. The stabilizing influence of 
the background temperature gradient is very marginal, since the core itself is close to being 
isothermal. Thus one may also expect to fold regions of very low density ratio Ro (also 
necessary for the development of the collective instability). We note, however, that the first 
He flash partially removes the electron degeneracy while the burning proceeds (see Figure 
[9]), so that one can only find a fingering region at relatively high Prandtl number for a very 
short time. Subsequent flashes take place in weakly degenerate conditions with Pr and r of 
the order of 10 -4 , which is too low for the collective instability to exist. As a result, the 
collective instability can at best only be active in a very short timeframe, and is therefore 
unlikely to play a significant role in stars undergoing core He flash. 


4.3. White dwarfs 


As for main sequence stars, fingering convection is not typically thought to occur spon¬ 
taneously in white dwarfs. However, the accretion of high-// material on their surface, from 
planetary debris or from winds from stellar companions, could trigger the double-diffusive 


instability, as recently discussed by Deal et al. (2013). 


Since electron degeneracy is high in a white dwarf except very close to the surface, we 
expect the Prandtl number and the diffusivity ratio to be relatively high as well, as in the 
cores of RGB stars. This is indeed the case, as shown in Figure [TOj For the purpose of 


illustration, we present a DA white dwarf model with effective temperature T e ff = 11150K, 
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kindly provided by G. Vauclair and S. Vauclair. This model was computed in the same 


manner as the DA white dwarf models of Deal et al. (2013), and its parameters were chosen 


to be representative of G29-38, the prototype DA white dwarf with a known debris disk 


and polluted atmosphere (Koester et al. 1997). Its mass is M* = 0.59M 0 , and its radius is 
A* = 0.0138A'.,. The Prandtl number and the diffusivity ratio are calculated as in Section 


4.2| (except that we use A = 12 and Z = 6 in the core of the star). For the purpose of clarity, 
we also show them as functions of the mass coordinate starting from the surface (i.e. with 
m = 0 corresponding to the surface, and m = M* corresponding to the core). Both ratios 
are quite small in the non-degenerate atmosphere (down to m/M* ~ 10 —6 ), then increase 
with depth to be of the order of a few times 1CP 3 in the rest of the star, suggesting that the 
collective instability might possibly be excited there. 

As with RGB stars, having a relatively high Prandtl number and diffusivity ratio is 
not enough to guarantee the presence of the collective instability. We must first deter¬ 
mine whether the fingering convection excited by material accretion may penetrate down to 
m/M* ~ 10 -6 and whether the density ratio at this level is indeed low. In order for the 
first condition to be satisfied, the mean molecular weight of the accreting material must be 
relatively large, at least larger than that of the non-degenerate envelope, and ideally larger 
than that of the core. This is not unlikely if its origin is from dust and planetary remnants. 
Whether a sufficiently low density ratio can also be created in the process remains to be 


determined with detailed calculations along the lines of those presented by Deal et al. (2013) 
for DA white dwarfs. This will be the subject of a follow-up investigation. 


5. Conclusion and observational prospects 

In this paper, we have shown through theory and numerical experiments that fingering 
(thermohaline) convection in stars can in principle drive large-scale gravity waves through 
a mean-field instability called the “collective instability”. Note that by large-scale we imply 
waves whose typical wavelengths are much larger than that of individual fingering modes, 
although from a global stellar point of view they remain quite small-scale (see below for 
more on this topic). We have performed a comprehensive exploration of parameter space 
to determine under which conditions the collective instability is excited, and concluded that 
it is limited to stellar regions in which the Prandtl number is relatively large (at least 
as large as 10~ 3 ). This only happens when electron degeneracy is high and dominates 
thermal conduction and viscous dissipation, such as in the cores of post main sequence stars, 
and white dwarfs. Compounded with the fact that the system must first and foremost 
be undergoing active fingering convection, this only leaves two possibilities for plausible 
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scenarios where the collective instability conld play a role: RGB stars undergoing core He 
flashes (albeit only for a very short time after the onset of the first flash) and white dwarfs 
undergoing accretion of high-metallicity material from planetary debris or companion stellar 
wind. 


The excitation of the collective instability in these examples could have two different 
consequences. First, it is a natural source of gravity waves - one could therefore wonder 
whether they may be directly observable through asteroseismology or not. Second, it causes 
a significant increase in the mixing rate (especially if convective layering occurs), which could 
in turn affect the evolution of the star, its internal structure and its surface abundances. 


As discussed in Section 4.2, the collective instability is only likely to be relevant in RGB 
stars in the first few tens of years after the onset of core He flash. The effect of enhanced 
mixing by fingering convection alone has (to our knowledge) never been investigated in detail, 
so we can only speculate on the added effects of mixing by collective modes. Both processes 
would likely act to erode compositional gradients established at the lower boundary of He 
flash convection regions. Whether this might have any observable signature, however short¬ 
lived, is unclear. Like all red giants, these stars exhibit “mixed modes” which behave like 


acoustic modes in the convective envelope and like gravity modes in the radiative core (Mosser 


et al. 2011). Their frequencies generally depend on the buoyancy frequency in the core, and 


are thus sensitive to composition gradients there. However, even without additional mixing, 
the magnitude of N^ (the component of the buoyancy frequency associated with //-gradients) 
established in He flash regions is much less than that at the outer boundary of the He core. 
The latter largely dominates the mode spectrum and overwhelms any asteroseismic signature 
coming from the He flash regions themselves (Hu et al. 2009). Since extra mixing can only 
reduce the small composition gradients resulting from He flashes, it would be difficult to 
deduce the presence of the collective instability from an observed mode spectrum. Finally, 
since //-inversions in the degenerate core only appear below He flash convection regions, the 
collective instability cannot modify surface abundances during this phase. 


The possibility of exciting gravity waves through the collective instability in white dwarfs 
is potentially more interesting on the other hand for two reasons. First, and as studied by 


Deal et al. (2013), fingering convection alone has a significant effect on the redistribution 


of high-// material accreted onto DA white dwarfs, and must be taken into account when 
estimating the accretion rate from observations of the surface chemical abundances. The 
additional mixing caused by the collective instability would go further in the same direction, 
and must therefore also be taken into account. Second, these waves could potentially be 
observable if they reach sufficiently large amplitudes and horizontal scales. For Pr ~ r ~ 
1CT 3 and reasonably small density ratios, the horizontal and vertical wavelengths of the 
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fastest-growing collective modes are of the order of 300d (see Figure [3]), which in dimensional 
terms correspond to about 30 meters. While this is too small to be observable, we note 
that larger-scale modes also grow, albeit at smaller rates (see Figure [l]). Furthermore, the 
nonlinear interactions between growing modes can also excite very large-scale ones even if the 
latter would otherwise intrinsically decay. In other words, it is not unlikely that global-scale 
waves could be present as well. The question of the saturation amplitude of the collective 
instability, however, is much more complicated, and remains to be addressed. Nevertheless, 
these findings raise an interesting prospect, namely that of finding polluted white dwarf 
pulsators outside of the standard instability strips (e.g. see 
review on white dwarf pulsators). 


Winget & Kepler 2008, for a 


Finally, we have confirmed the findings of Brown et al. (2013) concerning the sponta¬ 


neous formation of thermo-compositional layers as a possible nonlinear progression of the 
collective instability. Why these layers form and under which conditions one may expect to 
find them remains to be determined. However, this could mean that the ultimate signature 
of the collective instability may not be through the excited waves, but through the presence 
of thermo-compositional layering. 
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Fig. 1.— Largest value of the real part of A among the three possible solutions of the cubic 


(16), as a function of vertical wavenumber k (horizontal axis) and horizontal wavenumber L 


(vertical axis), for various values of the Prandtl number (with r = Pr) and reduced density 
ratio r = (-R 0 —1)/( 7_1 — 1)- Areas shown in white have negative growth rates. The “bulb” in 
each plot corresponds to the fingering modes, while the 
modes. 


‘leaves” are the collective instability 
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Fig. 2.— Properties of the fastest-growing collective mode (solid lines with circles), in 
comparison with the properties of the fastest-growing fingering mode (plain solid lines), as 
a function of the reduced density ratio r. Top left: real part of A max , the growth rate of the 
fastest growing collective mode, and of Af gm (solution of equations (19) and (20)), the growth 
rate of the fastest growing fingering mode. Top right: Oscillation frequency of the fastest- 
growing collective instability mode. Bottom left: horizontal wavenumber L max of the fastest 
growing collective mode. Bottom right: Vertical wavenumber fc max of the fastest-growing 
collective mode. 
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Fig. 3.— Total kinetic energy in the domain, and kinetic energy in selected collective modes. 
The modes are referenced by the number of wavelengths that fit in the computational domain, 
ffence the (1,0,1) mode has k 101 = (2-k/L x , 0, 2n/L z ), while the (2,0,2) mode has k 2 02 = 
(Air/L Xl 0, 4:ir/L z ). The theoretical growth rate of the fastest-growing fingering mode (see 
main text for detail) fits the observed exponential growth of the total kinetic energy at 
early times very well. The theoretical growth rate for the (1, 0,1) collective mode is Aioi = 
0.0094 ± 0.07 i, while A 2 o 2 = 0.02 ± 0.087 Both compare quite well to the corresponding 
collective mode growth observed in the simulation. By t = 300, the growth of both collective 
modes seem to have saturated, presumably due to their mutual nonlinear interactions. 
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Fig. 4.— Snapshots of the compositional held in our simulation at time t = 100 (left) and 
t = 350 (right). At the early time, the fingering held is fairly homogeneous, but at the later 
time, a collective mode with two wavelengths in each direction (the (2, 0, 2) mode) is clearly 
visible. 
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Fig. 5.— Compositional Nusselt number Nu M = Two layers form in the simulation, 

around t = 750. 
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Fig. 6.— Snapshot of the compositional held (left) and horizontally-averaged density profile 
(right) in our simulation at time t = 950. Two convective layers with nearly uniform density, 
separated by thin and strongly stratified interfaces, are clearly visible. Also plotted for 
comparison is the density profile at t — 350, showing that the gravity waves themselves do 
not affect the mean density profile much. 
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Fig. 7.— Left: Prandtl number v /kt as a function of radius in stars of various masses: from 
top to bottom following the arrow, M* = O.7M 0 , 0.9M 0 , 1.1M 0 , 1.3M 0 , 1.5M 0 , 1.7M 0 , 
2 Mq, 5 Mq, 1OM 0 and 30M 0 . Right: Diffusivity ratio k^/kt as a function of radius for the 
same stars. 
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Fig. 8.— Prandtl number (solid red line), diffusivity ratio (green dashed line) and degeneracy 
parameter C 2 (see text for detail, blue dotted line) for a 1 M 0 star just before the luminosity 
bump at t — 11.47Gyr. The left axis is for Pr and r, while the right axis is for C 2 . The 
shaded area mark the region of relatively high degeneracy (C 2 < 1) which is bounded from 
above by the H-burning shell. The values of Pr and r are least accurate around C 2 — 1. 
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Fig. 9.— Variation of the Prandtl number with mass coordinate as a function of time from 
the onset of the first core He flash (at t = 11.58Gyr) for a 1M 0 star. The entire core itself 
at that time spans the interval m G [0, O.45M 0 ]; only the region near the He-burning shell 
is shown. The line color is used as a proxy for time (from dark blue at the earliest time, to 
dark red at the latest time). Regions that are unstable to fingering convection are shown as 
a thick line in each profile. At the onset of the flash, and for a very brief period thereafter, it 
is possible to find fingering regions with high Prandtl number. However, the flash partially 
lifts the electron degeneracy, which causes the Prandtl number to decrease rapidly. 
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Fig. 10.— Prandtl number (solid red line), diffusivity ratio (green dashed line) and degener¬ 
acy parameter C 2 (see text for detail, blue dotted line) in the white dwarf model described 
in Section 43 In both panels, the left axis is for Pr and r, while the right axis is for C 2 . The 
shaded and white areas mark the region of relatively high (resp. low) degeneracy (C 2 < 1 f° r 
the shaded areas). The values of Pr and r are least accurate around C 2 = 1. Left: Profiles 
as a function of radius. Right: Profiles as a function of mass coordinate starting from the 
surface at m = 0. 









